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by 
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ABSTRACT 

A balloon borne measiirement of the cosmic ray electron spectrum 
from 10 to 200 GeV is reported in which two new techniques have been 
used to remove proton background contamination. First, the depth of 
the spectrometer was more than 25 radiation lengths , the equivalent of 
more than 2 mean free paths of material, enabling hadronically and 
electromagnetically induced cascades to be differentiated for a sub- 
set of the data. Second, electromagnetic cascade starting points were 
determined to within +0.1 radiation lengths based upon a calibration 
with electrons from 5.U to l8 GeV at the Stanford Linear Accelerator, 
greatly reducing the chances for a proton to simulate an electron. The 
resulting spectrum, when fitted with a power law, is quite steep, 

-3. 2+0.1, but the chi-square fit is marginal. A significantly better 
fit is achieved assuming a transition region model in which the source 
spectral index is 2.7 with a break occurring at about 50 GeV. 

^Portions of this paper will be submitted in partial fulfillment of 
the requirements for the Ph.D. degree from the University of Maryland. 
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Introduction 

The energy spectrum of primary cosmic electrons above 1 GeV has 
been given much attention in recent years because the electrons represent 
an important probe into the history and structure of the galaxy. They 
may provide information on the origin of all cosmic rays and on the back- 
ground of electromagnetic radiation. The electron spectrum is expected 
to exhibit a steepening at an energy related to the photon and magnetic 
energy density in the interstellar region where electrons are stored. 

While results up to several hundred GeV have been obtained in balloon 
flights, [Nishimura et al., 1969, Anand et al. , I968, 1969? Scheepmaker, 

1971j Meyer and Muller, 1971] serious discrepancies in the spectral shape 
and absolute intensity are apparent, and hence it is difficult to arrive 
at specific conclusions about the electron propagation. 

In the hope of clarifying the experimental situation, a large area 
detector was developed for use in a series of balloon flights. The ex- 
periment worked well on flights 1 and 3 which were flown from Alamogordo, 

New Mexico on April 28, 1969? at 7-9 g/cm^ residual atmosphere and on 
November 11, 1970 at 7.^ g/cm^. A total exposure factor of >=4000 M^-ster-sec 
was obtained. The experiment provided very detailed information on each 
event. Twenty three detectors were used to identify each event and sample 
the longitudinal electromagnetic cascade development. Each of these de- 
tectors was pulse height analyzed over a dynamic range of 10^ . A digital 
spark chamber provided trajectory data for each particle. Data on the 
primary cosmic ray proton and helium spectrum using this instrument has 
previously been published [Ryan et al., 1972]. 
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DESCRIPTION OF THE DETECTOR 

A diagram of the detector as flown on flight 3 is shown in figure 1 . 

It consists of three sections: a section for particle identification, an 
electromagnetic shower section, and a nuclear cascade section. The charge 
module is made up of four detectors; two l/ 4 " plastic scintillators, a 
1/2" acrylic plastic Cerenkov detector, and a 1/8" cesium iodide scint- 
illator. Each of these detectors is viewed by photomultipliers and their 
outputs are pulse height analyzed. Included in the charge module also 
is a digital spark chamber of the type developed by Ehrmann et al,, [1967]. 
This chamber allows particle trajectories to be determined so that geomet- 
rical corrections can be made to the pulse heights from the large area 
detectors. The spark chamber is also very useful in eliminating back- 
groimd particles entering the sides of the detector which satisfy the 
trigger criteria. 

The electromagnetic shower section of the experiment is composed of 
twelve modules, each of which consists of a I/8" thick tiingsten sheet 
and l/ 4 " plastic scintillator representing O.9 radiation lengths. Each 
module is viewed by two photomultiplier tubes whose outputs are summed 
and pulse height analyzed. 

Although the nuclear cascade section was primarily for measurement 
of the high energy proton spectrum, it was extremely useful for determining 
the interacting proton background correction to the electron intensity. 

The nuclear cascade section consists of seven iron modules, each one 
about 4.5 radiation lengths thick [0.5 nuclear mean free paths]. Again 
a sandwich structure is employed by alternating layers of plastic scintillators 
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and iron. Each iron module contains three plastic scintillators viewed 
by two photomultipliers tubes whose outputs are siommed and pulse height 
analyzed. 

TRIGGERING MODES OF THE EXPERIMENT 

In order to reduce backgroTind and the dead time in the experiment, 
the electron trigger criteria were set to demand the equivalent of a 7 
GeV electron in the tungsten modules. In addition a minimum pulse height 
was demanded in each of the 12 tungsten modules which required that 
the cascades develop very rapidly. This not only discriminated against 
low energy particles, but also eliminated triggers from high Energy pro- 
tons interacting deep in the tringsten. The only protons which triggered 
the experiment in the electron mode were those which interacted early 
enough in the tungsten stack to deposit a large energy in all the tungsten 
modules. 

Two other important modes were included in the experiments: these 

are the proton and calibration modes. In the proton trigger mode, a 
particle must deposit greater than Uo GeV in the nuclear cascade section 
of the experiment. Events were allowed to trigger in both the electron 
and proton modes simultaneously. These events were used to understand inter- 
acting proton background as discussed later. The calibration mode permits 
any charged particle which passed through both of the plastic scintillators 
in the charge mod\ile to trigger the experiment. This mode is used on 
the ground for calibration of individual detector elements on sea level 
muons. During the 1970 flight, every 32nd event opened a gate which per- 
mitted a calibration mode event to be recorded. In this way detector per- 
formance was monitored during the flight. 
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CALIBRATION OF THE EXPERIMENT 

While the instrument provides a great deal of information about each 
particle which triggers it, the transition effect [Piiikau, I 965 ] and the 
lack of any simple and reliable approximation to showers in heavy mater- 
ials make interpretation of the data difficult- It was realized early 
in the detector development that calibration of the instrument would be 
mandatory - 

The detector was calibrated on positrons from 5-4 to I 8 GeV at the 

Stanford Linear Accelerator, SLAG, in 1969 . The mean energy deposition 

for various energies and angles of incidence was normalized in terms of 

the energy deposition due to sea level muons which could penetrate the 

85 gm/ cm^ of the detector. All of these were then fitted to the expressions 

f (EQ,tQ,t)=exp (at+bt^) t^+0-5 (l) 

f eirp l-t+2 [(t-a) (y-0.56)]^/^ t>a+0.5 (2) 

where f (EQ,tQ,t) is the energy deposition that an electromagnetic cascade 

initiated by an electron of energy E^ will produce at depth t, if it 

started to cascade at depth to- The variables y and t are related to 

the incident energy and the depth by 

y=ln (|o), and t^ SD+lb 
^ COS0 

respectively where S is the detector module number, D is the depth in 
radiation lengths of one tungsten module, and 9 is the angle of incidence 
with respect to the normal, Ec,!),!^, and a are free parameters which are 
determined by simultaneously fitting the observed showers at eight dif- 
ferent energies to equations 1 and 2, [see Crannell et al., 1972 for a 
discussion of the calibration procedures!. The parameters a and b are 
uniquely determined by requiring f (Eo,tQ,t) and its first derivatives to 



be continuous at t=a+0.5« Equation (2) is after Heisenberg [19^6] and 
is based on approximation A of cascade theory. Equation (l) is simply 
used to carry the Heisenberg expression below the singularity at t=o,. 


Knowing Ec, D, and the energy deposition of an electron of unknown 
energy impinging on the instrument? the data from an event can be fitted 
to the electromagnetic cascade curves to determine an "apparent starting 
point" of the shower, t^, and the incident energy, E^. 


These parameters are determined by the minimization of 

■ y^= E [ f(Eo,to,ti) - fqi^ 
i=l 


(3) 


where cr^(EQ 5 tQ,tj^) is the variance of f (EQ,tQ,tj_) . At SLAG energies, 
we observed a^(Eo,tQ,tj_)i=4f (E^,tQ,tj_) except at small depths where the 
deviations were larger. \ also gives a measure of the likelihood of 
the incident particle being an electron and the data from SLAG provides 
a X distribution for electrons to compare against flight particles. 

Note that equation (3) represents a fit with 10 degrees of freedom so 
that Eq and to are well determined. 

The energy distribution observed by applying the procedure outlined 
above to a run with l8 GeV electrons from SLAG is shown in Figure 2. The 
FWHM of the distribution is ~17?^. The small peak below 40 GeV is due to 
two l8 GeV electrons incident on the experiment within the resolving time 
of the experiment [a frequent occurrence at SLAG, but not a problem in 
flight ] . 


DISGRIMINATION AGAINST PROTONS 


Since the electron component of cosmic rays is at best only a few 
percent of the proton component, the most serious problem in measuring 
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electrons is to reduce to a negligible contribution those interacting 
protons which may masquerade as electrons. Nuclear interactions which 
produce a tt° meson may produce a shower that is, in practice, indistinguish 
able from a shower produced by an electron. 

Our study of accelerator electrons and inflight high energy protons 
has revealed that the method of fitting the showers to the equation above 
provides a sensitive means of discrimination against the proton back- 
ground. In figure 3 a distribution of "apparent starting points" is 
presented for a sample of l6 GeV electrons from SLAG. This histogram 
is representative as no energy dependence was noted in the 5.4 to l8 GeV 
range. It should be noted that the "apparent starting point" distribution 
peaks at about 0.4l r. 1. before the tungsten stack. This agrees quite 
well with 0.29 radiation lengths of material that are calculated to be in 
the charge module preceding the tungsten stack and suggests that this 
distribution's width represents real fluctuations in the production of 
the first high energy cascade photon. 

In contrast to the electrons, a group of singly charged particles was 
selected from the 1970 flight which deposited >40 GeV in the iron and 
>5 GeV in the tungsten. The condition that the energy deposited in the 
iron exceed the energy deposited in the tungsten was also imposed to get 
a set of inflight interacting protons. Their "apparent starting point" 
distribution is plotted in figure 4. 

The distribution in figure 4 exhibits marked structure. The source 
of the structure is not completely understood. One quarter of the nuclear 
interactions should take place in the plastic scintillators which are 
invisible on this radiation length scale. The structure is probably the 
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result of the transport of Y-rays from tt° decay into subsequent tungsten 
slabs. However, the process of fitting electromagnetic shower curves to 
nuclear interactions may shift the "apparent starting point" .by a few 
tenths of a radiation length. In any case, the structure is real, is 
present for protons and especially for background events, and says that 
the shower fitting procedure is defining an "apparent starting point 
with an accuracy of approximately 0.1 radiation length. This has the 
effect of greatly reducing the amount of material in which a proton can 
interact to simulate an electron. 

ELECTRON SELECTION 

Each possible electron event is examined to be sure that it is caused 
by a singly charged particle as determined by the top plastic scintillator 
and the Cerenkov detector. Both of these detectors are 30 cm above the 
tungsten stack. If the particle is singly charged, the spark chamber data 
is used to extrapolate the path of the particle to the twelfth tungsten 
module. Only particles whose trajectories pass within the inner 80*)^ of the 
area of the 12th tungsten module are analyzed so a well defined geometry 
is used for flux calculations. 

Each event accepted is then fitted to the electromagnetic shower 
function to determine its energy, Eq, and apparent starting point, t^. 

The minimum \ is examined to see if its value is less than the 0.95 
probability point for electrons of energy Eq. If it is, and if the starting 
point is in the range for electrons observed at SLAC, the particle is 
accepted as an electron. 

RESIDUAL PROTON BACKGROUND AND OTHER CORRECTIONS 

Even though unambiguous identification of each incident particle 
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should be attainable with nuclear emulsions, emulsion techniques are re- 
sponsible for perhaps the highest [Anand et al., 1968 ] and the lowest 
[Marar et al., I 971 ] measured primary cosmic ray electron intensities. 

Other workers have used statistical separation of the protons and elec- 
trons. Scheepmaker [1971 ] has used the spectrum of interacting alpha 
particles in flight to estimate the interacting proton spectrum. Earl 
et al., [ 1972 ] and Meyer and MOller, [1971 1 have used data and extrapola- 
tions from accelerator pions and protons to determine the proton in- 
duced contamination. 

The \mique feature of this experiment is its total depth of four 
nuclear mean free paths of material. This depth allows secondaries from 
nuclear interactions to interact again and reveal themselves as proton 
induced showers. 

Although our rejection of such events is still a statistical one, 
we can use flight data itself to determine the proton background for 
this experiment. A subset of the apparent electrons [those which satisfied 
3-11 of the above criteria 1 is selected which have trajectories which pass 
through at least three of the seven iron modules. This subset of about 
2 ^ of the events is examined to determine the fraction which deposited 
significantly more energy in the iron than could be expected of an elec- 
tron of energy E. This leads to a 29+5^^ correction which is then applied 
to all the data. 

The final spectrum includes corrections to the data for energy loss 
due to bremsstrahlung in the residual atmosphere, atmospheric secondaries, 
spark chamber inefficiency and dead time due to paralysis of the experi- 
ment during event readout. The correction for re-entrant albedo electrons 
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is negligible at the energies measured. 

ENERGY DEPENDENT EFFECTS 

Extensive investigations for any energy dependent effect in selec- 
tion were also carried out. Energy dependence due to backscatter into 
the spark chamber or into the charge determination section was studied 
by examining the data from the spark chamber and the S2 scintillator 
below it, as a function of energy. No significant correlation was found 
between larger numbers of sparks and large pulse heights in the scintillator 
as a function of the estimated energy of the incident particle. No de- 
termination of any backward transport of gamma rays could be made, but 
the work of Ford [1972] indicates that back- scattering gamma rays greater 

than 2 MeV occur in less than 2^ of the electron showers at high energies. 

2 

A significant energy dependence was found in the y distribution, an 

2 

effect that was also observed in our positron data from SLAG. The y 
value below which 95^ of fhe SLAG positrons fell was found to increase 
like The spreading of the y^ distribution is qualitatively in 

agreement with our expectation since we observed that the standard devia- 
tions in the showers were larger than expected early in the showers. 

As energy rises there are more data points early in the shower [relative 

to shower maximum] and hence our underestimation of the deviations woixld 

2 

make the parameter y rise as a function of energy. 

To further check the data for energy dependent effects, the single 
particle events that failed the electron selection criteria were sorted 
into energy bins. Their spectral index was found to be -2. 7+0. 2, the 
same as the high energy proton spectrum [Ryan et al., 1972]. 

The final correction for proton contamination as determined by the 
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spectrometer is nearly energy independent, which is surprising because 
of the difference in spectra. Probably the proton rejection due to starting 
point location and cascade shape is better at high energies because the 
electromagnetic showers are less subject to fluctuations than their 
hadronic counterparts. This proves the background rejection ratio is 
imoroving at the same time the proton to electron ratio is going up. 

Vis-a-vis other experiments, we would like to point out that a pro- 
ton contamination of only 30 fo at low energies could become much larger 
at high energies since the electron and proton spectra are not the same. 
DISCUSSION 

While the intensity in figure 5 agrees with the measurements of 
many other workers at about 10 GeV, the spectrum is significantly steeper. 
Assuming the data can be represented by a single power law over the range 
10-200 GeV, the spectral index obtained is -3. 2+0.1. This spectral index 
is steeper than the often quoted value of ~2.7, [Bleeker et al., 19^8, 
Scheepmaker, 1971], although steep spectra have previously been reported 
[Earl et al., 1972, Nishimura et al,, 1969, Rubstov and Zatsepin, I 968 , 
Webber and Rockstroh, 19721. 

The spectrum shows no sharp "break" but does exhibit a consistent 
drop off in intensity. The Chi-square test for 8 degrees of freedom 
indicates that a power law represents a marginal fit to the data (al- 
most 3 standard deviations). 

The spectral index is expected to steepen by a full power if 
synchrotron and inverse Compton losses are dominant, based on the leakage 
lifetime approximation. This approximation has, however, been called 
into question by Jokipii and Meyer [1968] because it replaces the 
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diffusion term in the continuity equation by a characteristic loss term. 
Even the assumption of isotropic diffusion with a disk-shaped source 
region, produces a spectrum with two "breaks" each of one half a power 
at energies separated by a factor of 50 or more depending on the dimen- 
sions of the source region. Jokipii and Meyer predict a break in the 
1-10 GeV range, but they also point out that the assumption of isotropic 
diffusion is probably not realistic. 

Berkey and Shen [1969] have carried out a similar treatment using 
convection diffusion and find only one "break" occurs. However, their 
calculation does not include perpendicular diffusion nor the role of 
random magnetic field lines allowing escape from the disk. Using our 
data, the different models could certainly not be distinguished. A fit 
to a leakage lifetime model was thus attempted because of the poor fit to 
the power law. 

From Ramaty [l972], we find that the spectral break region can be 
represented in the following manner. If the source spectrum is of the form 

_p 

Q,(e)“'E'’ , then the equilibrium intensity I(e), would satisfy the expres- 
sion „ 1 /n o 

I(e)E^-Jq ^ (l-ox)^"^e"^dx (4) 

1 1 1 ^ 
with where W is the electromagnetic energy density in eY / car , and 

T is the mean electron lifetime in mega-years. The energy at which a=l 

is defined as the "break" energy. 

The transition region model fits the data at about the 1 standard 
deviation level and yields a value of 2.7 for F, the source Spectral index 
and a break energy of ~50 GeV. Wot too much significance should be attached 
to the particular model used or to the values obtained for the parameters; 
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the transition region model contained one more free parameter than the 

power law model, and the minimum in the chi-square was very broad. 

. 3 dir 

The two fits to the data are illustrated in figure 6 where E 
is plotted as a function of energy. [This representation is chosen to 
emphasize details of the data]. The best power law and the transition 
region model fits are indicated. The transition region model hreah 
corresponds to a life of about 6x10° years assuming W~1 eV/cm-^. 

Our data argues strongly against a spectrum as hard as -2. 7 around 
100 GeV, and is suggestive that the data cannot be represented by a single 
power law over the range 10-200 GeV. The theoretical interpretation of 
this result in terms of the leakage lifetime approximation is not clear. 

In view of the apparent complexity of the electron spectrum in the 
high energy region, new data above 100 GeV and more precise data below 
will be required to further our understanding. This is essential if 
we are to answer more questions about the cosmic electrons, their 
sources and their propagation through the galaxy. 
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FIGUEE CAPTIONS 


Figure 1 - Diagram of the ionization spectrometer showing the physical 
layout of the various component parts. 

Figure 2 - Estimated energy distribution for cascades initiated by l8 GeV 
electrons. 

Figure 3 - "Apparent Starting Point" distribution for electromagnetic 
cascades initiated by l6 GeV electrons. The origin of the 
abscissa references the top of the tungsten stack. 

Figure 4 - "Apparent Starting Point" distribution for singly charged 
cosmic ray primaries depositing ^0 GeV in iron stack and 
> 5 GeV in tungsten stack. 

Figure 5 - Differential spectrum of cosmic ray electrons (this experiment 
and other workers) 

Figure 6 - Differential spectrum of cosmic ray electrons at the top of 
the atmosphere multiplied by E^. The two successful flights 
are shown separately along with best fits to models. 



REFERENCES 


Anand, K, C., R. R. Daniel, and S. A. Stephens, Cosmic-ray Electron 
Spectrum above 50 BeV and its Implications for Cosmic-ray Con- 
finement, Phys. Rev. Lett,, 20, ?64, I 968 

Anand, K. C. , R. R. Daniel, and S. A. Stephens, The Cosmic-ray Electron 
Spectrum at Energies >100 GeV, Paper presented at the Xlth Int, 
Conf. on Cosmic Rays, Budapest, 19^9 

Berkey, G. B. and C. S, Shen, Origin and Propagation of Cosmic Ray 
Electrons, Phys. Rev., I 88 , 199^5 19^9 


Bleeker, J. A. M, , J. J. Burger, A, J. M, Deerenberg, A. Scheepmaker, 

B. N, Swanenburg, and Y. Tanaka, Energy Spectrum of Primary Elec- 
trons in the Range 0.5 - 30 GeV, Can. J. of Physics, S522, I 968 

Crannell, C. J, , R. A. Gearhart, F. A. Hagen, W. V. Jones, R. J. Kurz, 

J. F. Ormes, R. D. Price, R. F. Silverberg, and G. M. Simnett, 
Electron Calibration of a High-Energy Cosmic Ray Detector, Submitted 
to Nucl. Instr. and Methods, 1972 

Daniel, R. R. , and S. A. Stephen, Recent Observations on Cosmic Electrons 
and Their Consequences on Different Cosmic Ray Models and Related 
Astrophysical Quantities, Proc. Ind. Acad, Sci, 6 ^, 319> 19^7 

Earl, J, A., D. E. Neally, and T. A. Rygg, Balloon Measurements of the 
Energy Spectrum of Cosmic Electrons Between 1 and 25 GeV, J. Geo- 
phys. Res. 77, IO 87 , 1972 

Ehrmann, C, H, , C. E. Fichtel, D. A. Kniffen and R. W. Ross, A Magnetic 
Core Digitized Spark Chamber for Space Science Experiments, Nucl. 
Instr, and Methods, IO 9 , 19^7 

Faneslow, J. L. , R. C. Hartman, R. H. Hildebrand and P. Meyer, Charge 

Composition and Energy Spectrum of Primary Cosmic Ray Electrons, • 
Astrophys. J., I 58 , 771, I 969 


Ford, R., Private Communication, 1972 

Heisenberg, W. , Cosmic Radiation , Dover Publications, New York, pp, 11-25, 
1946 

Hovestadt, D., and P. Meyer, The Geomagnetic Cut-Off at Fort Churchill 
and the Primary Cosmic Ray Electron Spectrum from 10 MeV to 12 
GeV in I 968 , presented at the Xlth Int. Conf. on Cosmic Rays, 
Budapest, I 969 

Jokipii, J, R., and P, Meyer, Storage and Diffusion of Cosmic-Ray Elec- 
trons in the Galaxy, Phys. Rev. Letters, 20, 752, I 968 



REFERENCES (cont.) 

Marar, T. M. K,, P. S. Frier and C. J- Waddington, Intensity and 

Energy Spectrum of Energetic Cosmic Ray Electrops, J. Geophys. 

Res., 76, 1625, 1971 

Meyer, P. and D. MtCLler, Flux and Energy Spectrum of Cosmic Ray Elec- 
trons between 10 and 1000 GeV, presented at the Xllth Int. Conf. 
on Cosmic Rays, Hobart, Tasmania, 1971 

Nishimura, J,, E. Mikumo, I. Mito, K. Niu, I. Ohto, and I. Tiara, 

Observation of the Electron Component in Cosmic Rays with Emulsion 
Chamber, Presented at the Xlth Int. Conf. on Cosmic Rays, Budapest, 

1969 

Pinkau, K. , Errors in Electromagnetic Cascade Measurements due to the 
Transition Effect, Phys. Rev., 139 ? B1548, I965 


Ramaty, R., Private Communication, 1972 

Ryan, M. J., J. F. Ormes, and V. K. Balasubrahmanyan, Cosmic Ray Pro- 
ton and Helium Spectra Above 50 GeV., Phys. Rev. Letters, 28, 

985, 1972 

Rubstov, V. I., and V. I. Zatsepin, The Electron Energy Spectrum in the 
Flux of Galactic Cosmic Rays, Can. J. of Phys., 518, I968 

Scheepmaker, A. , Primary Cosmic Ray Electron Spectrimi between 5 and 
^ 300 OeV , unpublished thesis, University of Leiden, 1971 

Webber, W. R. and J. M Rockstroh, Cosmic Ray Electrons With E>lGeV-Some 
New Measurements and Interpretations, in press, J. Geophys. Res., 

1972 



-36.1 


t- 


1- 

5 


5 

Q. 


Q. 


■LUCITE LIGHT 
PIPE 



-108 


PMT 


P.M. 


P. M. 

nI 

P. M. 

si 

P. M. 

si 

A 

RM. 

M 

A 

P.M. 

n1 

y 

RM. 

*— — — S 


y/////////////////////////7?. 

y/////////////////////////7? 


y///////////////////////7Z. 


y///////////////////////ZZ 




y////////////////////////} 


y//////////7/////////Z77Z, 


//////////77////////////Z. 


’/////////////////////////. 


'//////////////////////TT?^ 


v//y///yyyy////////////z. 


v//////////////////////y^ 


y//////////////777y///77? 


y///////////////////////. 


'////////////////////////. 


y///////y/y///77//////77. 


v///////y/y//////////./7A 


v////////y////////y///77?. 

////////// / / / / J y / / / / ^ //-A 


•3mm TUNGSTEN 


■6.4mm PLASTIC 


■6.4mm PLASTIC 
SCINTILLATOR 


27 mm IRON 


IONIZATION SPECTROMETER 

th 



























zov. 



GY (GeV) 










n 


DEPTH IN RADIATION LENGTHS FROM 1st TUNGSTEN MODULE 



ELECTRONS/ (M'^-SEC-SR-GeV) 


□ FANSELOW et al ( 1969) 

▲ SLEEKER et al (1969) 
o HOVESTADT a MEYER (1969) 
O DANIEL et al ( 1967) 

V AN AND et al (1968) 

4 AN AND et al (1969) 

• SCHEEPMAKER (1971) 

■ NISHIMURA et al (1969) 

O MEYER a MULLER (1971) 

® PRESENT EXPERIMENT 



100 1000 
ENERGY (GeV) 




ENERGY IN GeV 


